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Recently Scholtz and Geyer proposed a very efficient method to compute metric op- 
erators for non-Hermitian Hamiltonians from Moyal products. We develop these ideas 
further and suggest to use a more symmetrical definition for the Moyal products, because 
they lead to simpler differential equations. In addition, we demonstrate how to use this 
approach to determine the Hermitian counterpart for a Pseudo-Hermitian Hamiltonian. 
We illustrate our suggestions with the explicitly solvable example of the — a; 4 -potential 
and the ubiquitous harmonic oscillator in a complex cubic potential. 
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1 Introduction 

Many non-Hermitian Hamiltonians H are known to possess real discrete spectra, 
e.g. ^12], which make them potential candidates for physical systems. Unlike as for 
Hermitian Hamiltonians h, the conventional inner products of the corresponding 
wavefunctions are usually indefinite and the central problem in this context is to 
construct meaningful inner products serving to formulate a consistent quantum 
theory. Besides from compensating wrong signs by hand jo], there are essentially 
three different, albeit in many cases equivalent, ways to achieve this: i) by employing 
bi-orthonormal eigenstates 0], ii) by constructing the so-called C-operator jS] |S] 
or iii) by restricting to pseudo-Hermitian Hamiltonians H and constructing their 
Hermitian counterparts h related to H by a similarity transformation -q — [5] 

h = rjHr]- 1 = /i f W = tfHrf 2 . (1) 

See also [§] for an earlier discussion of these issues. In many respects the last 
possibility is the most direct and straightforward approach on which we will almost 
exclusively concentrate here. The natural starting point is usually a given non- 
Hermitian Hamiltonian H for which one needs to construct rj. Subsequently one 
may formulate all relevant physical questions for the non-Hermitian system in terms 
of the conventional formulation of quantum mechanics associated to the Hermitian 
system. Thus a key task in this approach is to find rj 2 and rj. Unfortunately this is 
only possible in some exceptional cases in an exact manner and otherwise one has 
to resort to perturbation theory. Building on earlier work |dT ll(J|. we found in |11| 
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the closed expressions 

rli ri±ii 

l 2l Inl |\nr 1 2 I 2n— 1,, 

h = ho + J2 9 4 i7 2 L 4 2 "^o), F = Ao - £ 7^T)T C ^" )(/l0) ' (2) 

n— 1 ^ ^ n— 1 ^ ' 

which are related according to Q perturbatively. Here [a;] denotes the integer part 
of a number x, the non-Hermitian Hamiltonian is assumed to be of the form 

H = ho + igh\, (3) 

ho = fig, hi = h\ and Cg"'(x) denotes the n-fold commutator of the operator q 
with some operator x. In case one has the condition Cq +1 \ho) = for some finite 
integer i the expressions (0 are exact and otherwise they are just correct up to the 
stated order in the coupling constant g. The E n are Euler's numbers and the K2n-i 
may be constructed from them One should note, however, that, in practice, the 
above-stated procedure may lead to rather cumbersome relations involving commu- 
tators. This fact poses a major problem for determining closed formulae for specific 
isospectral pairs of Hamiltonians, or for carrying out perturbative computations up 
to higher orders. 

The main purpose of this paper is to propose a practical scheme for overcoming 
this difficulty. We elaborate further on a recent proposal by Scholtz and Geyer ^21 
to solve by means of Moyal products instead of computing commutators. The 
central idea is to exploit isomorphic relations between commutator relations and 
real valued functions multiplied by Moyal products, which correspond to differential 
equations. We shall demonstrate that this approach is rather practical and allows 
to compute pairs of isospectral Hamiltonians h = h) and H ^ H ' . 



2 Similarity transformations from Moyal Products 
2.1 Generalities on Moyal products 

Moyal products are applied in a wide field of research, such as non-commutative(nc) 
M-theory 15 , nc-string theory JHjj nc-integrable field theories |17l 1181 ITill I2()j . 
etc. The key idea is to transfer the noncommutative nature of some operators to 
real valued functions. Technically one may set up such an isomorphism in various 
different ways. In the present context of studying non-Hermitian Hamiltonians such 
possibilities have been exploited in ^2). The authors defined the Moyal product of 
real valued functions depending on the variables x and p as 

f(x,p) * g(x,p) = f(x,p)e^*g(x,p) = -d s J{x,p)d s p g{x,p), (4) 

The classical, more widespread and symmetrical definition is, see e.g. ^] |^ E2] 

f(x,p)*g(x,p) = f(x,p)eH d * d *-*d^g(x,p) (5) 
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or exchange the roles of x and p in J2J. In order to achieve a proper isomorphism 
between the operator expressions and those computed with Moyal products one 
requires different types of prescriptions to translate real valued functions into oper- 
ator valued expressions. Products computed with definition Q must be viewed as 
ordered products in which all p-operators are moved to the left of all i-operators. 
In expressions computed with definition JSJ on the other hand one should replace 
each monomial p m x n or x n p m by the totally symmetric polynomial S mn in the m 
operators p and n operators x 



Here ir indicates the sum over the entire permutation group. The simplest example 
to illustrate this is 



We use here the standard canonical commutation relation [x,p] = i and throughout 
the paper we employ atomic units K = e = m e = cot = 1. We observe that the *- 
product yields the correct operator expression upon replacing x — > x, p — > p. The 
★-product on the other hand corresponds to the correct commutation relations for 
px — > = (px + xp)/2. The defining relation @ is slightly more complicated 
than Q , but it is more symmetrical and leads therefore to cancellations of various 
terms as one can easily convince oneself. Loosely speaking, it has the advantage 
that it incorporates already more of the noncommutative nature of x and p from 
the very beginning. As we shall see below, this is the reason why it leads to simpler 
differential equations for the quantities we wish to determine. These statements 
are also supported by comparing the equations resulting from (UJ) and (@J) for all 
examples calculated up to now with this method ^2 1 ) ■ 

2.2 Hermitian counterparts from Pseudo-Hermitian Hamiltonians 

Let us now briefly discuss in general how we proceed to construct h from a given 
H by first solving Q for r/ 2 and r/. The explicit knowledge of r\ is vital, since once 
it is known one may control the entire quantum mechanical formalism, such as 
inner products, observables, time evolution, etc. However, in general one can not 
compute rj exactly and has to rely on perturbative methods [HJ 1101 II 1 1 in which 
one has to solve the commutator relations occurring in J2J) order by order. This is 
a very cumbersome procedure and up-to-now it has only been carried out for few 
cases to lowest order. Here we present a simple and more efficient scheme which 
leads to the exact determination of rj by employing Moyal products. We build on 

1 ) For instance for the example H = p 2 + igx 3 , definition yields 



2gx i r) 2 (x,p) + 3igx 2 d p rj 2 (x,p) — 'igxd 2 rj 2 (x, p) — igdpr) 2 (x, p) + 2pd x ri 2 (x , p) + id 2 rj 2 (x,p) = 0, 




(6) 



[x ,p \ =Aipx~2 = x *p —p *x =4ipx — 2 = x -kp —p *x' 



2 



Aipx. (7) 



whereas J^J gives the simpler form &gx 3 r) 2 (x, p) — 3gxd 2 rj 2 (x,p) + 4,pd x r) 2 (x,p) = 0. 
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suggestions of jTJ], but as a starting point we use instead of the definition (@J the 
definition ijjjj and write the second equation in Q isomorphically as 

* rj 2 = rj 2 * ff. (8) 

With definition © this will then lead to a differential equation in rj 2 , whose order 
depends on the degree of x and p in H. In most cases this equation can not be 
solved exactly and under these circumstances we assume a perturbative expansion 
for rj 2 in the coupling constant g, which was introduced in @ 

V 2 (x,p) = V°° g n c n {x,p). (9) 

* — 'n— 

Combining JHJ and then leads to a differential equation, which involves the 
functions c„ in a recursive manner and can therefore be solved order by order. Ul- 
timately we aim at exact expressions, which yield rj 2 to all orders in perturbation 
theory Having then solved various differential equations, we naturally expect 
some ambiguities in the general solutions, which mirror the possibility of different 
boundary conditions. We would like to stress that this is not a drawback, which is 
only present when using Moyal products, but it is rather a reflection of a general 
feature of perturbation theory. The same kind of ambiguity occurs in the perturba- 
tive approach based on commutators. In that context one may only fix the operator 
corresponding to the c n (x,p) up to any operator which commutes with the Hcrmi- 
tian part of H, that is Hq. This means that, in the expressions are insensitive 
to any replacement q — > q + q with [q, ho] — 0. A further type of ambiguity, which 
is always present irrespective of an exact or perturbative treatment, is a multipli- 
cation of rj 2 by operators which commute with H, i.e. we could re-define r/ 2 — > rj 2 Q 
for any Q, which satisfies [Q, H] = 0. 

We may fix these ambiguities by invoking a further property of 77. As discussed 
in assuming a dependence on the coupling constant g of the form 

v(-9)=v(gr\ h(g) = h(-g) and rf(g) = H{-g), (10) 
will guarantee the pseudo-Hermiticity relations . Therefore we require next 

V 2 (9)*ri 2 (-9) = 1, (11) 

which may be solved systematically order by order when we already know the 
expansion up to the ambiguities. Note that (|II|I is automatically satisfied if we 
take rj 2 = e q and q = X^^Li ff 2 ™ -1 92n-i as assumed in many cases on the grounds 
of PT-invariance [5]. Having determined r/ 2 the transformation i] is subsequently 
computed easily order by order from 

V*V = V -V n 9 n c n (x,p), (12) 
when assuming a further power series expansion 

Eoo 
g n q n (x,p). (13) 
n— 
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Finally, having obtained an explicit expression for 77, we may compute the Hermitian 
counterpart to H from (JJJ as 

h*r) = rj*H. (14) 

According to the above arguments, i.e. the second equation in (JTUJ, we should find 
an expression of the general form 

h(x,p) = Y]™ g 2k h 2k (x,p). (15) 

Let us now illustrate the above with two explicit examples, by starting with one for 
which the differential equation may be solved exactly and thereafter complicating 
it to a case which requires perturbation theory. 



3 The non-Hermitian -x 4 potential 

As a straightforward example we consider the non-Hermitian Hamiltonian 

H(xJ)=p 2 - P - + a(x 2 -l) +z.g (^lp--2 ax y (16) 

which results from H = —cP/dz 2 — ez 4 when using z = —2i\J\ + ix as transforma- 
tion, a = 16e and the introduction of the coupling constant g to separate off the 
non-Hermitian part |2.'{j . The exact similarity transformation r\ for this Hamilto- 
nian was recently constructed by Jones and Mateo |23 using perturbation theory in 
terms of commutators in the spirit of equation (J2J) . See also j2SJ for further reason- 
ing on how the Hamiltonian (|16|) can be used to make sense of the — ez 4 -potential 
despite its unappealing property of being unbounded from below. 

In order to illustrate the method let us see how to determine rj by using Moyal 
products. First we notice that we wish to treat H(x,p) as a real valued function 
and we therefore have to replace the anti-commutator with the appropriate Moyal 
products. When using definition © we have to replace {x,p 2 } by x*p 2 +p 2 *x = 
2xp 2 and the differential equation JHJ for the Hamiltonian Ijl6|) reads 

= Agp 2 xr] 2 {x,p) — 8gxarj 2 (x,p) — 4xad p r] 2 (x,p) (17) 
-d x r] 2 (x,p) + 4pd x r) 2 (x,p) + 2gpd p d x rj 2 {x 1 p) - gxd 2 x J] 2 [x,p). 

As a comparison we also present the differential equation resulting from JSJ when 
using the *-product instead of the *-product. When converting the operator valued 
Hamiltonians into a function, we have to pay attention to the fact that (/ * g)* ^ 
f* * g*. Thus in H we replace i{x,p 2 } by i(x *p 2 +p 2 * x) — 2ip 2 x —2p, whereas in 
we substitute —i{x,p} by — i(x*p 2 +p 2 *x) — ~2ip 2 x +2p, which is of course not 
the same as converting first the anti-commutator with a subsequent conjugation. 
The resulting differential equation reads 

= 4gp 2 xrj 2 (x,p) ~ 8gxari 2 (x,p) — Axadpri 2 (x,p) (18) 
+i2gp 2 d p r] 2 (x,p) - i2ad 2 r/ 2 (x,p) - iAgad p r] 2 (x,p) + i4gpr) 2 (x,p) 
-(1 + 2g - Ap - Apxig)d x r 1 2 (x,p) + (2i - 2gx)d 2 x rr 2 {x,p). 
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Obviously, equation 1(18(1 is more complicated than 1(17(1 , which illustrates our asser- 
tion that the ★-product is more advantageous than the ^-product. Ultimately they 
should lead, however, to the same result. Indeed, each line in (|17|) as well as (|18|) 
vanishes separately for 

ry 2 (x,p) = e^- 2 ^, (19) 

which, when compensating a slight difference in convention, is precisely the same 
expression as found in |2rij . To find rj{x,p) from r] 2 (x,p) by means of (|12|l is trivial 

9 P 3 

in this case as we just have to take the square root, i.e. rj{x,p) = e~s~^~~ 9P . Using 
1(14(1 thereafter we find the Hermitian counterpart of H 

h{x,p,g)=p*-?-+a(x 2 -l) +g Ap2 ~ 2a)2 . (20) 
2 4a 

When setting the artificially introduced parameter g to 1, we recover precisely the 
expression found in 

4 

h(x,p,g=l) = ^--^+x 2 a. (21) 
4a 2 

In also the interesting massive case H = —d 2 /dz 2 + m 2 z 2 — ez A has been 
discussed. The transformation from the z-variable to the rr-variable will add in ((1611 
a term — m 2 (l + Aigx). The resulting similarity transformation is then rj(x,p) ~ 
e TTr-9P( 1 + 2m2 / a ) anc i the Hermitian counterpart reads 

h(x,p,g)=p 2 -^ + a(x 2 -l) -4m 2 +g 2 ^ >- 22 

I 4a 

which for g — > 1 reduces precisely to the expression reported in |2.'{j . 

Let us next embark on an example for which the differential equation l|17Jl can 
not be solved exactly. 



4 Harmonic oscillator perturbed by a complex cubic potential 

The prototype example for the study of non-Hermitian Hamiltonian systems is 
the harmonic oscillator perturbed with a complex cubic potential 

tf = Y + y + ^ 3 - (23) 

It was the discovery yQ that this Hamiltonian possesses a positive real discrete 
spectrum, which led to the current interest in this subject. This Hamiltonian is 
obviously non-Hermitian, but PT- invariant and pseudo- Hermitian pi IH1 171 ITH] . The 
latter property means that the relations Q hold. For the Hamiltonian 1)23(1 the 
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conjugation relation JHJ translates into the differential equation 2 ) 

= 4gx 3 f] 2 (x,p) - 3gxd 2 r, 2 (x,p) - 2xd pV 2 (x,p) + 2pd x i 1 2 (x , p) . (24) 

Using the expansion 112|) we find the recursive equation 

= 4x 3 c„_i(a;,p) - 3xd 2 c n -i(x,p) - 2xd p c n (x,p) + 2pd x c n {x,p). (25) 

Notice the occurrence of the aforementioned ambiguities. It is easy to see that we 
may add to c n {x,p) any arbitrary function ([(p 2 + x 2 )/2]. This type of ambiguity 
just reflects the fact that in the perturbative formulation in terms of operators we 
may add to q any function of h a = (p 2 + x 2 )/2. We write c n (x,p) — c n (x,p) + 
C[(p 2 + x 2 )/2] and determine first the functions c n (x,p) which also obey 
Taking as the initial condition co(x,p) = 1, we find order by order 



- / s 4 P 3 „ 2 - / s , 2 2p 2 x i 8 a; 6 
ci{x,p) = — + 2px , c 2 {x 1 p)=6x — , 

112p 5 64p 9 56p 3 x 2 32p 7 x 2 A . 56p'V 140p 3 a; 6 

c-i(x.p) = 1 V\vpx 

iy ,P> 15 81 3 9 y 9 27 

I6px s 

2 2 128gV no 4 1856pV , 32p 8 x 4 5392p 2 x 6 
C4{x,p) = IVlp x h 98a; 



3 15 81 45 

416p 6 .x 6 1768a; 8 230p 4 x 8 176p 2 a; 10 160a; 12 
243 45" + 81 + 81 + 243 ' 



c 5 (x,p) 



oo lB 3 . 2752p 7 19456p 
P 5 405 


11 : 2048p 15 
' 3645 


- 336px 2 + 


9632p 5 a; 2 
5 


107008p 9 a; 2 1024p 13 a; 2 


+ 1928p 3 x 4 - 


26752p 7 a; 4 


3328p n a; 


405 243 


45 


243 


8332px 6 272896p 5 x 6 


18304p 9 x 6 


154132p 3 a; 8 


760/a; 8 


15 405 


729 


405 


+ 27 


11488pa; 10 t 23524p 5 a; 10 


5536» 3 x 12 


320pa; 14 




135 1215 


H 

729 


+ 243 ' 





This can be expanded effortlessly to higher orders, but we will not report these 
functions here. Next we demand the dependence on the coupling constant g to be 
of the form (|ll(l in order to fix the ambiguities. Since we may add to c a {x,p) any 
arbitrary function £q [ho — (p 2 + x 2 )/2] , due to the recursive equation (|25|) this 



2 ) For comparison definition Jl) yields 

= 4g x 3 r) 2 (x, p) — 2 x dpfj 2 (x, p) + 6 i g x 2 d p r) 2 (x, p) — id 2 rj 2 (x, p) 
-6gxd 2 r/ 2 (x,p) - 2 i g d^rj 2 (x , p) + 2pd x r) 2 (x,p) + i d 2 r?(x, p). 
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function will produce descendents at higher level in c n for n > a. Thus we have the 
general form 

c n (x,p) = c n (x,p) + C*°M- ( 26 ) 

We can now use this function and achieve that is satisfied order by order. 
Explicitly, 

*(..P) = ^ + 

C2(a;,p) = c 2 {x,p) + Co 2) [^o] = c 2 (x,p) + -fr 3 , - 4/i . 

As we note, indeed the additional function just depends on ho. In a similar fashion 
we can compute the higher order functions 

248» 5 32n 9 6Ap 3 x 2 16p 7 x 2 , 8n 5 x 4 4p 3 x 6 

4 832/ 32p 12 2 2368 P 6 64p 10 2 

ct^p) = 152 P ___ + __ + (56p -^^ + — )x 

v 3 9 7 v ^ 9 ; 3 ' 

r ^ 109^3^2144/ 4672p n 128p 15 ,10864/ 
c 5 (x,p) = -1024p H — — h + (— 168p 



5 405 3645 v r 15 



li 



20416/ 64p 13 2 712/ 3488/ 64p 
405 + ^43~ )X + ( ~ 45~ + ~M~* X 

, 32/ N „ ,28/ 8/ N o 4/x 10 
-(28p + 48/ - — ^)x 6 - (— ^- tt)x ' 



27 7 v 3 9 7 15 ' 

„„„,, 2 98336/ 340256/° 6016/ 4 256/ 8 , „„„ 

.ufiifin 4 70592p 8 35392/ 2 256/ 6 , 2 2 70288/ 

+6160p - + - — )x - (216p - —g- 

26816/° 64/ 4 . 4 . 1520/ 3232/ 128/ 2 , R 

+ T)5 W^ + (340 + -3 45- + "24T )a; 

, 2 328/ 16/°, o ,20/ 16/ N 10 4// 2 

We observed that always ^ 2 ™ +1 ^ [/j ] = o. Making next the ansatz (|13|l . we compute 
the q n {x,p) by reading off the corresponding powers in the expression 112(1. We find 
order by order 

2/ 2 , 2 2p 6 x 2 2p A x 2 p 2 x A 

qi(x,P) = ~ir + P x > Q2{x,p) = -p- + — + — 



3 ' 7 9 2 3 2 

79/ 4/ 23/x 2 2pV 13px 4 /x 4 p 3 x 6 



(Ja fa;, pi = 6p 

wv r 15 81 3 9 4 3 6 

, s 67/ 148/ 2/ 2 ,37b 2 442/ 4/° % 9 

fl4 I.p = 1 h ( 1 JSC 

y v 2 45 243 v 2 45 81 1 
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,61 29p 4 
■(-5- + -T- 



P \ 4 



7P 2 



_ p_ 

9 

405 ' 3645 



)X' 



-997p 3 443p 7 434p 



li 



1952p y a;^ 2p vi x z 1199p J x 4 



405 
149px 6 
+ — 
1677p 2 
2 

16811» 4 
+ ^2- 
1361p 10 
+ 405 
943p 2 
~32 



243 

- 6p 5 x e -r 

103649p 6 
90 

8987p 8 



p x 



12 

43p 3 a; 8 



24 ! 
355p:r 2 
" — 4 
361pV 
-15- 4 



5^,2 



4681p 5 x 

w 

2p n x 4 



+ 



7 8 

p X 



p X 



81 

10 



27 24 36 

84373p 10 286/j 14 



1350 
1717p 12 



1215 

4pi6 



120 ' 
4p 18 
32805 
2503p 2 



1215 
1861 



+( 



45 
1731,8 + 



16 



72 



108' 



3645 
1625p 4 

12 
29p 4 



359p 8 2p : 



1131 

4 

43021/ 



180 



90 



243 



)x 6 



48 



180 ; 720 



Finally, having obtained an explicit expression for the rj, we may compute the 
Hermitian counterpart to H from 114|) . We find order by order 



h (x,p) 

h 4 (x,p) 
h 6 (x,p) 



El 

2 



x 

Y 



1 - 3i 4 
hz{x,p) = -- + 35 2 ,2 + — 



27p 2 + 2p b 



15i 2 



128 - 984p 4 - 72p 8 
27t 8 

+2885 2i6 + ^—, 



— 36S , 4 i 2 - 
14645, 



515 



2.4 



2 

7685 6i2 



7i 6 



(27) 



84 5; 4 + 66054 



where we have already converted from real valued functions to operators. Notice 
that h(x,p) is of the form (|15f) and as we expect all terms h%k—i — 0, such that 
indeed the second relation in i|10|) holds. Up to hi(x,p) our results agree precisely 
with |24j. Higher orders do not seem to appear in the literature. 



5 Conclusions 



We have shown that Moyal products can be used as a very powerful tool to com- 
pute the similarity transformation 77 relating a pseudo- Hermitian Hamiltonian H to 
its Hermitian counterpart h. As concrete examples, we have applied such products 
to a situation for which 77 can be calculated exactly, namely the non-Hermitian — x 4 
potential, and to a case for which it can only be computed perturbatively, that is 
the harmonic oscillator with a complex cubic perturbation. In the latter case, we 
profited considerably from the fact that the Moyal products involve only real- valued 
functions, instead of commutators, and lead to differential equations which can be 
solved recursively. Relation was crucial to fix the ambiguities, which result 
from operators which commute with the unperturbed Hamiltonian h . 
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In a more general context, this is a far less cumbersome procedure than those 
involving commutators, and allows one to evaluate isospectral pseudo-Hermitian- 
Hermitian pairs perturbatively up to a very high order. Even more, Moyal products 
can be employed to obtain closed formula for the similarity transformation rj to all 
orders in perturbation theory. Such formulae have been obtained for several specific 
isospectral pairs, and will be discussed in a subsequent publication |13| . 

Despite the practical use of this approach, there are clear limitations, such as 
for instance when the potential is of a nonpolynomial nature. In that case the sum 
in the definition for the Moyal product does not terminate and one has to deal 
with differential equations of infinite order. Hence, there is a clear need for further 
alternative methods, such as recently proposed in |25j . 

Acknowledgments. We are grateful to Hugh Jones for useful comments. 
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